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Abstract 

We present the next-to-leading order QCD corrected total cross sections and the distributions 
of the transverse momenta of the final anti-bottom-quark, top-quark and charged Higgs-boson 
for the processes of pp/pp thH^ + X in the minimal supersymmetric standard model(MSSM) 
at the Tevatron and the LHC. We find that the NLO QCD corrections significantly modify the 
leading-order distributions of the transverse momenta of final particles(p^, p\. and p^ ), and the 
total NLO QCD corrections reduce the dependence of the cross section on the renormalization 
and factorization scales, especially the NLO QCD corrected cross sections at the LHC are nearly 
independent of these scales. Our results show that the relative correction is obviously related to 
rriH- and tan /?, and the total NLO QCD relative corrections can be beyond —50% at the Tevatron 
and approach —40% at the LHC in our chosen parameter space. 



PACS: 12.60.Jv, 14.80. Cp, 14.65.Fy 
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1 Introduction 



One of the major objectives of future high-energy experiments is to search for scalar Higgs bosons and 
investigate the symmetry breaking mechanism of the electroweak interactions. In the standard model 
(SM) rr, one doublet of complex scalar fields is introduced to spontaneously break the symmetry, 
leading to a single neutral Higgs boson H^. But there exists the problem of the quadratically divergent 
contributions to the corrections to the Higgs boson mass. That is the so-called naturalness problem. 
The super symmetric (SUSY) extensions of the SM provide a possibility to solve this problem. In 
supersymmetric models, the quadratic divergences of the Higgs boson mass can be cancelled by loop 
diagrams involving the supersymmetric partners of the SM particles exactly. The most attractive 
and simplest supersymmetric extension of the SM is the minimal supersymmetric standard model 
(MSSM)[21 O]- In this model, there are two Higgs doublets Hi and H2 to give masses to up- and 
down-type fermions. The Higgs sector consists of three neutral Higgs bosons, one CP-odd particle 
{A^), two CP-even particles {h^ and H^), and a pair of charged Higgs bosons {H^). 

However, these Higgs bosons haven't been directly explored experimentally until now. If the 
charged Higgs boson is lighter than the top quark, it would probably be found at the upgraded 
Tevatron or at the future LHC through the t bH^ decay channeljl]. Otherwise, if the charged 
Higgs boson is heavier than the top quark, there are three major channels to search the charged 
Higgs boson: (l)Charged Higgs boson pair production El El E]; (2)Associated production of a 
charged Higgs boson with a W boson^l^J; (3) Associated production of a charged Higgs boson with 
a top quark^2|. The decay of the charged Higgs boson has two major channels: H~ ^6|13j. and 

The associated production of a charged Higgs boson with a top quark (pp/pp tH^ +X) seems to 
be the most promising channel [T^. Corresponding to ignoring or observing the final state anti-bottom 
quark experimentally, the cross section can be inclusive (ignore final state anti-bottom quark) or ex- 
clusive(observe final state anti-bottom quark). The inclusive and exclusive leading order subprocesses 
can be written as: 

gb^tH-, (1.1) 

gg^tbH-, qq^tbH~, (1.2) 

respectively. The next-to-leading order(NLO) total cross section of the inclusive process has been 
studied in the SM QCDj^ and supersymmetric 0CD[T7j. 

In this paper, we calculate the total cross section of exclusive processes for the associated production 
of the charged Higgs boson with top quark and anti-bottom quaik{pp/pp — > tbH~ + X) in the MSSM 
at hadron colliders including the NLO QCD corrections. In section 2, we present the calculations of 
the leading order cross sections to pp/pp — > tbH~ + X in the MSSM. In section 3, we present the 
calculations of the NLO QCD corrections to pp/pp tbH~ in the MSSM. The numerical results 
and discussions are presented in section 4. Finally, a short summary is given. 
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2 The Leading Order Cross Sections 



The exclusive tbH^ production mechanism at the parton level contributing to the hadronic processes 
pp/pp tbH~ + X, involves qq {q = u, d) annihilation and gluon-gluon fusion channels. The subpro- 
cess via the qq annihilation is written as 

q{pi)q{p2) ^ tip3)biP4)H-ip5), (2.1) 

And the subprocess via gluon-gluon fusion is denoted as 

9iPi)9iP2) ^ tip3)b{p^)H-{p5), (2.2) 

In above notations of these two channels, we use pi, p2 and p^, p^, p^ to represent the four-momenta 
of the incoming partons and the outgoing particles respectively, and write them in brackets. The 
Feynman diagrams at leading order(LO) for the two subprocesses (|2.1j) and (|2.2j) are plotted in Fig.l 
and Fig. 2, separately. In Fig.l both figures are gluon exchanging s-channel diagrams with a charged 
Higgs boson being radiated from anti-top-quark and bottom-quark, respectively. Fig.2(a-b), Fig.2(c-e) 
and Fig.2(f-h) belong to the s-channel, t-channel and u-channel diagram groups, respectively. 




Fig. 2. The tree- level Feynman diagrams for gg — > tth subprocess. 
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For the ()2.1() subprocess qq — > thH , the corresponding LO and NLO amplitudes can be expressed 
in the form as|19j: 

^LO,NLO ~ ^'^'^^'lO,NLO (2-'^) 

where C^^ is the only color factor involved in the LO amplitude of the subprocess qq thH" , which 
can be written as: 

= » A^ (2.4) 

The first 3x3 SU{3) Gell-Mann matrix A'^ arises from qq color state, the second Gell-Mann matrix 
A'^ arises from tb color state. 

Similarly, the LO amplitude of the (|2.2|) subprocess gg — > tbH~ can be expressed as: 

Mfg = Qcf + Cf + Cf^?^ Mf^ + Qcf - Cf + Cf^'^ Mf , (2.5) 

with 

Mfa = + ^M|^ (2.7) 

Mf = M^-^Mf, (2.8) 

where Cnin = 1,2) are the color indexes of incoming gluons, /"^'^ and d"'^'^ are the SU{3) antisymmetric 
and symmetric structure constants respectively, matrixes 1 and A'^ arise from tb color state. Mf^, 
M^^ and Mu^ are the amplitude parts corresponding to the s-, t- and u-channel amplitude groups, 
respectively. Taking the trace of color factors, we get: 

on 

rr((:7f ^Cf ) = cf%fc with cf = 24, cf = 48, cf = — . (2.9) 

The squared amplitude for the LO can be written as: 

|Mfg|2 = — (|Mf |2 + |Mf |2) - — . 2i?e(Mf ^ • Mf ). (2.10) 
3 3 

Then the LO cross section for the subprocesses qq, gg tbH~ can be obtained by using the following 
formula: 



"LO 



I d^sEl^ir I'' (2-11) 
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where d$3 is the three-particle phase space element. The summation is taken over the spins and colors 
of initial and final states, and the bar over the summation recalls averaging over the spins and colors 
of initial partons. The LO total cross section of pp/pp tbH~ + X can be expressed as: 

a^o{AB{pp,pp) ^ tlH- +X) = 
(gg) ^ ^ 

^ 1 I A / dxAdxB[Gi/AixA, IJ')Gj/BixB, IJ')^'LoixA,XB, IJ') + iA^ B)], (2.12) 

ij=(uu),(aa) 

where xa and xb are defined as 

XA = -^,XB = -^, (2.13) 

where Pa and Pg are the four-momenta of the corresponding protons or antiprotons. A and B represent 
the incoming colliding hadrons(proton/antiproton). o^^loiij = uu,dd,gg) is the LO parton-level total 
cross section for incoming i and j partons. Gi/A{B)^^ are the LO parton distribution functions (PDF) 
with parton in a proton/antiproton. 



3 NLO QCD Corrections in the MSSM 

In the calculation of the NLO QCD corrections in the framework of the MSSM, we adopt the 't 
Hooft-Feynman gauge, and use dimensional regularization(DR) method in D = 4 — 2e dimensions to 
isolate the ultraviolet (UV), infrared(IR) and collinear singularities. Renormalization and factorization 
are performed in the modified minimal substraction(M S*) scheme. The NLO QCD corrections can be 
divided into two parts: the virtual corrections from one-loop diagrams, and the real gluon/light-quark 
emission corrections. 



3.1 Virtual Corrections 

The virtual NLO QCD corrections in the MSSM come from self-energy, vertex, box and pentagon dia- 
grams. We plot some QCD one-loop pentagon diagrams in Fig.3 for demonstration. These corrections 
in the MSSM can be sorted into two parts. One is the so-called SM-like QCD correction part coming 
from the diagrams with gluon/quark loops, another is SUSY-QCD correction part arising from virtual 
gluino/squark exchange contributions. The amplitude for the virtual SM-like NLO QCD correction 
part contains both ultraviolet (UV) and soft/collinear infrared (IR) singularities, while the amplitude 
corresponding to the NLO SUSY-QCD diagrams contains only UV singularities. 

In order to remove the UV divergences, we need to renormalize the strong coupling constant, the 
wave functions of the relevant fields and the masses of bottom-, top-quark. In our calculation we 
introduce the following counterterms. 

mt — > mt + Smt, mb^mb + SiJib 
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tL 
bL 



1 + -6Zl ] tL, tR 



1 + -6Zji ] tn 



1 + -6Zl j bL, bn ^ [1 + -SZ'j, j bn 
l + ^6Zl]qL, qn ^ (l + Uzl\qR 



(1 + -5Zg)G^, 



9s ^ gs + Sgs 



(3.1) 



where gs denotes the strong couphng constant, t, b, q, and denote the fields of top-, bottom-, 
up(down)-quark and gluon. The wave functions of the relevant fields, top quark mass in propagators 
and in the Yukawa couplings are renormalized in the on-shell(OS) scheme. For the renomalizition of 
the strong coupling constant gs, we adopt the MS scheme at renormalization scale fir, except that 
the divergences associated with the top quark loop and colored SUSY particle loops are subtracted 
at zero momentum (20) . Since we define the counterterm of the strong coupling constant consisting 
of SM-like QCD term and SUSY QCD teTTa{6gs = 5gi^^^ '^'^^^ -|- 5g^f^^^^), these two terms can be 
obtained as 



{SM-like) 



in 



r){SM—like) 

^ - + -ln 

2 e 3 



2 



(3.2) 



{SQCD) 



4tt 



-i \ In ^ 

2 e 3 //2 



=1,2 



+ 



where we have used the definitions: 



(SM-like) 
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-N 



If 



E 

U =u,c,t 



(SQCD) 
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In 



i=l,2 

+ E 

D=d,s,b 



m 



12 



In- 



//2 



-N- 



,(3.3) 



(3.4) 



The number of colors = 3, the number of light flavors nif = 5 and 1/e = l/ec/y — 7_e + ln(47r). 
The summation is taken over the indexes of squarks and generations. Since the MS scheme violates 
supersymmetry, it is necessary that the qqg Yukawa coupling gg, which should be the same with the 
qqg gauge coupling gg in the supersymmetry, takes a finite shift at one- loop order as shown in Eq. (|3.5|1 



9s = gs 



1 + 



Svr V3 



-N-Cf 



(3.5) 



with Cp = 4/3. In our numerical calculation we take this shift between gs and gs into account. 

In the Yukawa couplings, we use the MS mass of the bottom quark, rnb{fir), to absorb the large 
logarithms which arise from the renormalization of bottom quark mass(21). but keep the bottom quark 
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pole mass everywhere else. The bottom quark mass in propagators is renormalized by adopting the 
on-shell(OS) scheme. The expressions of the MS mass of the bottom quark rni,{^j.) corresponding to 
1-loop and 2-loop renormalization groups are given by: 



I— loop 
rnbiHr)2-loop 



rrib 



rrib 



as{mb) 
Ois{mb) 



co/bo 



CO /bo 



, Co , , . asilir) - as{mb) 

L + — [Ci — Oi) 

bo 



TT 



3 



TT 



where 



bo 
bi 



1 /^ll.. 
— —N nif 

47r V 3 3 ^ 



1 51N 



2 

—1 
3 

19nif 



Co 



1 

IT 



1 



Cl 



(lOliV - Wnif) 



2iT nN-2nif ' 72tt 
The renormalization of the bottom quark mass in Yukawa couplings is defined as: 



'-rub 



(3.6) 

(3.7) 

(3.8) 
(3.9) 

(3.10) 



where the counterterm for the SM-like QCD S'^^^ is calculated in MS scheme, while SUSY-QCD 
counterterm S^Q*^^ ig calculated in on-shell(OS) scheme. 

Since there are significant corrections to tbH~ production for large values of tan (3, we absorb these 
corrections in the Yukawa couplings j22j. The resumed thH~ Yukawa coupling can be expressed as: 



9tbH- 



V2f 
V I 



mt cotf3 Pr + mb{fir) tan/3 



1 



tan^ /3 



1 + 



(3.11) 



where |23| 



Ab 
Anib 
Ai 

I(a, b, c) 



Arub 

- rUg 12 tanP I {7711^,1711^,1711) 

-I ^"^9 I{ml,7nl,ml), 

a b c 
ab log — + bc log — \- ca log — 
b c a 

(a — b){b — c)(c — a) 



(3.12) 



When we use the resumed tbH Yukawa coupling to express the tree-level cross sections, we have 
to add a finite renormalization of the bottom quark mass in tbH~ Yukawa coupling to avoid double 
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counting in the NLO QCD cross section. '24V 



nib 



Af = ^ ^ (1 + ^^^) mg /i tan /3I(m?^,m|,m?). 

In the calculations of one-loop diagrams we adopt the definitions of one-loop integral functions 
as in Ref.,26^. The Feynman diagrams and the relevant amplitudes are generated by using FeynArts 
3|27j. and the Feynman amplitudes are subsequently reduced by FormCalc32[2S\ . The phase space 
integration is implemented by using Monte Carlo technique. For the IR-finite integral functions, the 
numerical calculations are implemented by using developed LoovTools^29\. 

The numerical calculations of the IR-infinite integral functions are implemented by using the 
methods described in Ref. j3Uj . The method is given by: 

rp(N)D _ rp(N)D I 1 _ I 

fj,i...fj,P — -'^ fj,l...fj.p\stng -r ^j^,„^P ^ ^i...^p\smg, {O.IO) 

where T^^]^p and T'^}^p\sing are the A^-point integral functions and their complete mass-singular 
parts in D = 4 — 2e dimensions, respectively. t'^}^^p and T^f^.^'.^p | ^ are in 4 dimensions. t'^}'^^p 
can be calculated by using mass rGnornictlization sclieme. Tji-^\^p\giyig can be calculated by: 



Tl^}.^,p{VO, ■ ■ ■ ,PN-l,mo, . . . ,niN-l)\sing = 
N-1 N-1 

^ ^ AnkC^^...^p{pn,Pn+l,Pk,mn,nin+3,mk)- (3.14) 
n=0 fe=o 

where C^i...;ip is the corresponding 3-point integral functions. C^^,,,^p in 4 dimensions can be calcu- 
lated by using mass renormalization scheme. The explicit expressions of the L^i...;jp in D = 4 — 2e 
dimensions can be found in Ref. |31j. 

The 0{as) QCD virtual corrections of the cross sections in the MSSM to the subprocesses qq, gg — > 
tth^ can be expressed as 



~ (99.99) 



I d$3E2Re (-^It^'^-^iS?) > (3.15) 



where M^rel -^irel the Born amplitudes for qq,gg tthP subprocesses, and -A/i^i^nai ^^'^ 
M^yiJtuai ^'■^ renormalized amplitudes of all the NLO QCD Feynman diagrams involving virtual 
gluon/quark and gluino/squark for qq annihilation and gg fusion processes, respectively. 

Then the amrtuai is UV finite, but still has IR divergence. Its soft IR divergence part can be can- 
celled by adding with the soft real gluon emission corrections, and the remaining collinear divergences 
are absorbed into the parton distribution functions. 
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3.2 Real Parton Emission Corrections 



The 0{as) corrections due to real parton emission give the origin of IR singularities. These singularities 
can be either of soft or collinear nature and can be conveniently isolated by slicing the phase space 
into different regions defined with suitable cutoffs, a method which has a general name of phase space 
slicing(PSS) 32 . In our calculation we consider the real parton emission subprocesses listed below for 
a consistent and complete mass factorization: 



q{Pi) + Q{P2) t{p3) +b{p4) + H'ip^) + g{pQ), 

g{Pi)+9{P2) t{p3.) +b{p4) + H~{p5) + g{p(i), 

{q,q){pi)+g{p2) ^ tfe) + + ^~(p5) + (g,g)(P6), {q = u,d). (3.I6) 

Using the method of phase space slicing method [H^. we introduce an arbitrary small soft cutoff 5s 
to separate the 2 — > 4 phase space into two regions, according to whether the energy of the emitted 
gluon is soft, i.e. < 6s\'^/2, or hard, i.e. Eq > 5sVS/2. In the real gluon emission processes 
qq — > tbH~g and gg tbH~g, there are both soft and collinear IR singularities. While in the real 
light-quark emission processes (in (|3.16|) ). there are only collinear IR singularities, but no soft ID 
singularity. The cross sections for the real gluon emission subprocesses listed in 1)3. 16|) can be divided 
into two parts to isolate the soft IR singularities: 



areaiiqq, gg ^ tbH g) = asojtiqq, gg ^ tbH g) 

+ ahardiqq,gg ^ tbH'g), (3.17) 

where agoft is obtained by integrating over the soft region of the emitted gluon phase space, and 
contains all the soft IR singularities. Furthermore, we decompose d'hard of the four subprocesses (in 
(|3.16|) ) with real gluon/light-quark emission, into a sum of hard-collinear (MC) and hard-non-collinear 
{HC) terms to isolate the remaining collinear singularities from a^ard^ by introducing another cutoff 
6c named collinear cutoff, i.e.. 



^hard{qq,gg,qg,q9 ^ tbH ig,g,q,q)) = aHciqq,g9,q9,qg ^ tbH ig,g,q,q)) 

+ ^Hc((l^^Sg,qg,q9 ^ tbH'{g,g,q,q)). (3.18) 
In the HC regions of the phase space the following collinear conditions are satisfied: 

^ n: < or < (5c, (3.19) 

and at the same time the emitted parton remains hard, anc contains the collinear divergences. The 
analytical expressions of the cross sections in the soft and HC region, asoft and anc, can be obtained 
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by performing the phase space integration in d-dimension. In the HC region, is finite and can be 
evaluated by using standard Monte Carlo techniques in 4-dimensions (M] . The cross sections of the 
real gluon emission process pp/pp — > tbH~ g + X, (Tgoft-, (^HC and o"-^, depend on the two arbitrary 
parameters, 5s and 6c- However, after mass factorization the total NLO QCD corrected cross section 
(^reai{pp/pp ^ tbH~ g + X) is independent on these two arbitrary cutoffs. We shall explicitly discuss 
that in Sec. 4. This constitutes an important check of our calculation. In the next two subsections, we 
will discuss in detail about the soft and hard-collinear gluon emission. 



3.2.1 Soft Gluon Emission 

For the real gluon emission subprocesses 

q{Pi)q{P2) t{p-i)b{pi)H~ {p^)g{pQ), g{pi)g{p2) t{ps)b{p4)H~ {p5)g{pe) (3.20) 

the soft region of the phase space is defined by 

< £^6 < SsVI/2. (3.21) 

The gluon bremsstrahlung cross sections of both qq and gg collision channels can be written in the 
following form|19j: 



^soft = '5"LO ® ^ XI ■ ^j)9^AP^^Pj)^ 



where Tj are the color operators jEllSEllS3, 9ij are the soft integrals defined as: 



9^j{Pi,Pj) 



(2vr/u) 



2e 



271" JEe<SsVS/2 Eq 



'^{PiPj) P| 

{PiP6){PjP6) {PiPiif {PjP&f 



(3.22) 



(3.23) 



The similar expressions of the soft integrals for qq,gg — > ttH^ in the SM can be found in Ref jl8l I19j . 

Using the definitions of color operators, we get the expressions of (Jsoft for qq annihilation and gg 
fusion channels, respectively. 



QQ 
soft 



as 
2^ 



17 1 

^{912 + 534) - ^(513 + 524) - :j(5i4 + 923) 
bod 



^QQ 



(3.24) 



gg 

soft 



^ I d^3E [(^^1 + 161^3) |Mf 1^ + (^Z)2 + 161)4) |Mf |2 



+ 



+ 161)3) 2i?e(Mf ^ • Mf ) 



(3.25) 
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where Mf^ and have been expressed in Eas. H2.7|) and (|2.8j) respectively, and the notations, Di, 
Z)2) and 1)4, used in above equation are defined as 

Di = 9gi2 + 9^13 + 9924 - 534, 
D2 = 9gi2 + 9^23 + 9gi4 - 934, 

D3 = 6(512 - 514 - 523 + 934), 

D4 = 6(512 - 513 - 524 + 534), (3.26) 

3.2.2 Collinear Parton Emission from the Initial Parton 

9. Hard Gluon Emission Subprocess qq{gg) tbH~g 

Let the hard gluon be emitted collinear to one of the incoming partons, i.e., 

-TT-?^ < '^c, or < 5c. (3.27) 

In this region, one of the initial state partons, i{i = q,q,g), is considered to split into a hard parton 
i' and a collinear gluon, i.e., i — > i'g, with pi' = zpi and pq = {1 — z)pi. The matrix element squared 
for qg{gg) tbH~g factorizes into the Born matrix element squared and the Altarelli-Parisi splitting 
function, and is expressed as: 

Y,\MHc{^j - tbH~g)\' (47ra,/i2^)^|M,o(i'j ^ tbR- g)\^ l' -2Pu'{z,e) \ ^g) 



where 



z 1 — z 

+ + z{l-z) 



1 - z 
1 + z^ 



Pggiz) = 0, 



P,,{z) = CF\^j--j, P^g(z) = -CF{l-z). (3.29) 

Using the approximation pi — p^ zpi{i = 1,2), the element in the four body collinear phase space 
region can be written as|18j 

d^4\coii = _ ^^ zdzdiiQ[-{l - z)iiQ]-'e (^ii^g^l^ - 5,6^ , (i = l,2), (3.30) 

where s' = 2pii ■ pj. Note that the four body phase space should be evaluated at a squared parton- 
parton energy of zs. Therefore, after integration over the collinear gluon degrees of freedom, we 
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obtain 



O'HC 



Us r(l — e) / AiTfi: 



27rr(l - 2e) V s' 





[(l-^)n 


/ dz 


/o 


2z 



e)aLo{i'j thH ) + {i ^ j) 



(3.31) 



In order to factorize the collinear singularity into the parton distribution function, we introduce a 
scale dependent parton distribution function using the MS convention: 

as r(l -e) / iTTfilY (I 



Gq/A{x,^if) = Gq/Aix) 



^ 27rr(l -2e) \ fi} J 

as r(l - e) /fTT^yi f^-^' 
27r r(l -2e)[ fi} J \ J, 

Us r(l - e) ( 47r/i^ Y / 1 



dz 



Af{q - qg) 



Pqq{z,e)Gg/A{x/z), 



(3.32) 



+ 



^ 27rr(l - 2e) I fij 

as r(l - e) / ATTfij 
27rr(l-2e) fij 



-) ^r(5-55)) 



^^^'^ dz 

—Pgg{z,e)Gg/A{x/z), (^ = p,pp.33) 



where 



Af{q^qg) = Ci.(21n,5, + 3/2), = 4/3, 

^r(5^59) = 27Vln5, + (lliV-2n//)/6. 



(3.34) 



By using above expressions, the NLO QCD correction parts of the total cross sections contributed by 
qq annihilation and gg fusion subprocesses in the initial state collinear phase space region are obtained 

as 
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HC 



^qq 



as r(l — e) / Anfif. 



27rr(l - 2e) V s 



,2\ <^ 



[Gq/AixA, fJ'f)Gg/BixB, fJ-f) + Gq/AixA, fJ'f)Gq/BixB, fJ'f) 



d,d 



Al'^{a ag) 



+ Al%a ag) 



Gq/A(.XA, fJ-f)Gg/B{xB, fJ-f) + {A^ B)\ dxAdxBi3.35) 



gg 
'Jhc 



-gg 
'Jlo 



as r(l — e) / 47r/i^ 



+ 



27rr(l - 2e) V « 

^^^^^ + K{g^gg) 



\ {Gg/A{^ 



2 y,^yXA,^J.f)Gg/B{xB,^^f) + Gg/A{xA, ^J'f)Gg/BixB, fj-f) 

Gg/A{xA, fJ-f)Gg/B{xB, f^f) + {A^ B)\ dxAdxB, (3.36) 
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where 



ASC 
^0 



Ga/A,B{x,l^f) 



dy 



Ga/A,B{x/y,iif)Paa{y), (« = u,u,d,d,g), 



with 



I (1 - y)'^ S \ _ - 

-Paa(y) = -Paa(y) In -| J ~ Kaiv)' {a = u,u, d, d, g). 



(3.37) 



(3.38) 



2. Hard Light-quark Emission Subprocesses {q, q)g tbH^ + {q.. q) 

The method in the calculation of the hard light-quark emission subprocesses {q, q)g tbH~ + (q, q) 
is similar to that for hard gluon emission subprocesses. In the collinear region, the initial state parton 
i{i = u, d, u, d, g) is considered to split into a hard parton i' and a collinear light-quark, i — > i'q, with 
Pi' = zpi and pq = (1 — z)pi. Let the hard light-quark be emitted collinear to one of the incoming 
partons, the collinear region is then defined as: 



2pi • P6 
EqVI 



or 



2p2 • P6 



< Sc. 



The collinear singularity of o-^^al '^^^ ^® written as: 



CThC 



Ug r(l — e) / Ann: 



,2\ « 



dz 



27rr(l-26) \ s' J \y e)^^^[Jo 2 J 

■P,9i^, ^)^IU99 - tlH-) + Pg,{z, e)dfo{qq ^ thH-)] . 



with 



(3.39) 



(3.40) 



Pii'{z,e) = Pu'{z) + ePl,,{z), 



1 /l + {l-zf 



P'gg{z) = -z{l-z), 
iV^ - 1 

2N \ z r ^99(^)= 
Using the MS scheme, the scale dependent distribution function can be written as: 



a, T(l-e) 



27rr(l-2e) \^ a*/ 



^ dz 



(3.41) 



PiAz)Gi,A{x/z). (3.42) 



14 



And we can get the expression for the initial state coUinear contribution at 0{as)^ order: 



i=q,q 



&fQ{xA,XB,l^) 

dz 



s (l-zfSc 



Pig (z) In 



JxA Z ^' Z ' 



d-fo{xA,XB,l^) 



3.3 Total NLO Cross Section 



s (l-zfSc 



Pgi{z)ln 



(3.43) 



The final result for the 0{as) NLO QCD corrected cross section part of qq annihilation subprocesses 
can be written as: 



'-'NLO 



J dxAdxBGqiA{xA, n)Gq/BixB, jj) ^loi^A, Xb, fJ-) + ^Ittuali^^, Xb,^) + ^HftixA, XB,fJ') 

+ {A^ B)] + a%% + J dxAdxB [Gq/A{xA, l^)Gq/B{xB, n)a^ixA,XB,l^) + (A ^ B)] . (3.44) 
And the NLO QCD corrected cross section part for gg fusion subprocess has the expression as: 

'^NLO ^ \ j dxAdxBGg/A{xA, lj)GgiB{XB,lj)[^fo^XA, XB,n)+ ^yfrtuali^^^ XB,lj)+ a^g^f^ixA, XB, n) 

+ {A^ B)] + ajfc + \j dxAdxB [GgjAixA, fi)Gg/BixB, n)al^ixA, xb,IJ') + {A ^ B)] (3.45) 
The cross section of {q, q)g — thH" + (g, q) {q = u, d) can be written as: 

^NLO = ^fc+m J dxAdXB[Gi/A{xA,l^)Gg/B{xB,l^)^^{xA,XB,l^) + (A^ B)]. (3.46) 



i=q,q 



with the hard-non-collinear partonic cross section given by 



HC 



Y}M{ij^tbH-g{q,q))\^d^i. 



(3.47) 



where (i$4 denotes the four-particle phase space element. 

Finally, the NLO QCD corrected total cross sections for pp/pp —>■ tbH~ + X can be obtain by 
using the formula: 

^11 _L _L ^19 (3.48) 



, gg I qg 

O-NLO — (^NLO + '^NLO + '^NLO- 



In our calculation, we have checked the cancellations of the UV and IR divergence analytically, and 
the final results are both UV- and IR-finite. 
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4 Numerical Results and Discussions 



In our numerical calculation, we adopt the CTEQ6M[3Sj parton distribution functions and the 2-loop 
evolution of as(/i) to evaluate the hadronic NLO QCD corrected cross sections with a^^^ {Mz) = 
0.118, while for the hadronic LO cross sections we use the CTEQ6L1 parton distribution functions 
and the one-loop evolution of as(/^)- We take the SM parameters as: aew(-^z)~^ = 127.918, niy/ = 
80.423 GeV ^ mz = 91.18 GeV, mt = 175 GeV, rrif, = 4.62 GeV. As a numerical demonstration, in 
this work we refer to the Snowmass point SPSlb for the relevant MSSM parameters |38j. if there is no 
other statement. The MSSM parameters in this benchmark are given by: 

tan/3 = 30, /i = 495.6 GeV, At = -729.3 GeV, 
Ah = -987.4 GeV, m-g = 916.1 GeV, rug^ = 762.5 GeV, 

= 780.3 GeV, mi^ = 670.7 GeV. (4.1) 

The charged Higgs boson mass is taken as free parameter. The calculations are carrying out at the 
upgraded Tevatron with the pp colliding energy y/s = 2 TeV and the LHC with pp colliding energy 
= 14 TeV. 

As the check of the correctness of our numerical calculation, we plot Fig. 4 to show the dependence 
of cr^fa^ on 6c and 5s at the LHC. Fig.4(a) shows the o"^f„^ as the functions of 5c by taking 5s = 10~^, 
mu- = 310 GeV and 5c varying from 10"^ to 10~^. Fig.4(b) presents the curves of cr^f^^ versus 5s 
with 5c = 10~^, rufj- = 310 GeV and 5s running from 10~^ to 10~^. Both figures show that our 
results of (y^eai^dd ~^ tbH~ + g) is independent of cutoff 5c and 5s. Actually, our calculation proves 
also that the real gluon emission cross section for qq annihilation channel is independent on cutoff 5c 
and 5s, and the cross sections for qg and qg fusion channels(as shown in (|3.16|l ) are independent on 5c 
too. 

In following calculation, we require the final anti-bottom quark to have a transverse momentum(p^) 
being larger than 20 GeV and a pseudorapidity |r/b| < 2 for the Tevatron and less than \r]h\ < 2.5 
for the LHC, unless other statement is given. We present the dependence of the cross section on 
the renormalization/factorization scale Q/Qo (where we take Q = fir = fJ-f for simplicity and define 
Qo = {mt + mfj- +mh)/2 ) in Fig.5(al-a2) and Fig.5(bl-b2) taking m^^- = 175 GeV at the Tevatron 
and rufj- = 250 GeV at the LHC, separately. In Fig.5(al) and Fig.5(bl) we depict the curves 
for total cross sections at leading-order aio and (Jnlo including NLO corrections for the processes 
pp/pp tbH~ -|- X at the Tevatron and the LHC. They show that the (Tnlo in the MSSM is less 
dependent on the normalization/factorization scale Q/Qo than a^o both at the Tevatron and the 
LHC, especially in Fig.5(bl) the curve for unlo at the LHC is nearly independent of Q/Qq. We 
can see from Fig.5(al) that the NLO QCD corrections suppress the LO cross section of the process 
pp tbH~ -|- X at the Tevatron. Fig.5(a2) shows that at the Tevatron the dominant contribution to 
(^NLO is from the subprocess qq thH~ , while Fig.5(b2) demonstrates that the main contributions 
are coming from subprocess gg — > thH~ at the LHC. From Fig.5(al) we can see that if Q goes down to 
a very low value, i.e., Q « Qq, the curve for unlo tends to have a negative value. The reason is that 
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large logarithmic corrections spoil the convergence of perturbation theory in the proton-antiproton 
colliding energy at the Tevatron. From Fig.6(al-a2) and (bl-b2), we conclude that the dependence of 
the NLO QCD corrected cross section ajVLO on the scale Q is significantly reduced comparing with 
(TlO: especially at the LHC, there the a^LO is very stable in a large range of Q. 

In the following figures, we fix the value of the renormalization/factorization scale being Q = Qo = 
IJ,r = Hf- In Fig.6(al) and (bl) we depict the LO and total NLO QCD corrected cross sections, cr^o 
and (TNLO-, in the MSSM as the functions of rriH- at the Tevatron and the LHC, respectively. The 
corresponding relative NLO QCD corrections 5{= ) versus rxiji- are plotted in Fig.6(a2) 

and (b2), separately. From these figures, we can see that the cross sections aNLO and cr^o decrease 
rapidly as ra^- varies in the range from 175GeV to 30QGeV at the Tevatron and from llbGeV to 
hbQGeV at the LHC. We can read from Fig.6(al) that when rrijj- increases from 175 GeV to 300 GeV , 
the total NLO QCD corrected cross section aNLO in the MSSM decreases from 0.66 fh to 0.02 fh at 
the Tevatron. From Fig.6(a2) wc can read when mjj- increases from 175 GeV to 550 GeV, (Jnlo 
decreases from 160 fh to 3 fb at the LHC. The absolute value of the relative correction 6 in Fig.6(a2) 
decreases slightly form 39% to 32% at the Tevatron. While the absolute value of the relative correction 
5 in Fig.6(b2) increases rapidly with the increment of the charged Higgs boson mass at the LHC, it 
reaches the value of 44% when m^- = 550 GeV . 

In Fig.7(a-b), we present the total LO and NLO QCD corrected cross sections a^o and unlo 
in the MSSM as the functions of the transverse momentum of anti-bottom quark cut p'^^y^^, taking 
niH- = 175 GeV at the Tevatron and m^- = 250 GeV at the LHC, respectively. In these two figures, 
we can see that the cross sections ctnlo and a^o decrease as Px cut go^s up from 10 GeV to 25 GeV. 
We can read from Fig.7(a-b) that when Pj^^ut increases from 10 GeV to 25 GeV , the total NLO QCD 
corrected cross section gnlo in the MSSM decreases roughly from 0.87 fb{55 fb) to 0.56 fb{37 fb) 
for the Tevatron(LHC). 

Fig.8(al) and Fig.8(bl) show the LO and total NLO QCD corrected cross sections aio and aNLO 
in the MSSM as the functions of tan/3 with rrijj- = 175 GeV at the Tevatron and nijj- = 250 GeV at 
the LHC, respectively. When tan /3 has a small value, the NLO QCD corrected cross section decreases 
from 0.45 /6(30 fb) to 0.15 /6(127 fb) at the Tevatron(LHC), as tan/? goes up from 4 to 12. The 
NLO QCD corrected cross sections unlo increases from 0.15 fh{\2 fb) to 2 fh (132 fb) for the 
Tevatron(LHC), as tan/3 varies from 12 to 50. The corresponding relative corrections 5 are plotted 
in Fig.8(a2) and Fig.8(b2). The 5 can be beyond —50% at the Tevatron, and approach —40% at the 
LHC at the position of tan/3 ~ 15. 

In Fig. 9, we depict the distributions of the transverse momenta of the final states(py, and p^ ) 
with rriij- = 250 GeV at the LHC and mjj- = 175 GeV at the Tevatron, separately. In Fig.9(al-a2, 
bl-b2 ,cl-c2), we show the distributions of the differential cross sections daLO,NLo/dPT^ daLO,NLo/dpj- 
and daLO,NLo/dp^ at the LHC and the Tevatron, respectively. These figures demonstrate that the 
NLO QCD corrections significantly modify the leading-order distributions of the differential cross 
sections daLo/dp\-, da^o/dpip and daLo/dp^ at hadron colliders. We find that in the low pj, region 
the differential cross sections daLO,NLo/dpj^ can be very large, especially at the LHC. 
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5 Summary 



In this paper we calculated the NLO QCD corrections to the processes pp/pp —>■ tbH~ + X in the 
MSSM at the Tevatron and the LHC. We investigated the contributions of the NLO QCD corrections to 
the total cross sections and the distributions of the transverse momenta of final particles (anti-bottom- 
quark, t-quark and charged Higgs-boson), and found the NLO QCD corrections significantly modify the 
corresponding LO differential cross sections. We analyzed the dependence of the NLO QCD corrected 
cross sections and the corresponding relative corrections on the renormalization/f actor ization scale Q, 
charged Higgs-boson mass rrifj-, transverse momentum cut of anti-bottom quark Pj^cMt' ^^'^ tan/3, 
respectively. Our numerical results show that the theoretical NLO QCD corrections in the MSSM 
reduce the dependence of the cross section on the factorization and normalization scales, especially 
the NLO QCD corrected total cross sections at the LHC are nearly independent of these scales, and 
the relative correction is obviously related to mjj- and tan/3 at both the Tevatron and the LHC. We 
find the total NLO QCD relative corrections 5 can be beyond —50% at the Tevatron and approach 
-40% at the LHC. 
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Fig. 4. The cutoff dependence of cr^^^i with Mjj- = 310 GeV at the Snowmass point SPSlb at the 
LHC. Fig.4(a) shows the Sc dependence with Sg = 10~^, and Fig.4(b) shows the 6s dependence with 
6c = 10-^ 
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Fig. 5. The total cross sections a^o and unlo for the processes pp/pp tbH~ + X as the functions 
of the renormahzation/f actor ization scale Q with rriH- = 175 GeV, \r]i,\ < 2 at the Tevatron and 
nifj- = 250 GeV, \r]i,\ < 2.5 at the LHC, are shown in Fig.5(al) and Fig.5(bl) respectively. The 
contribution parts to the total cross sections UNLoiPP/pP ~^ tbH~ + X) from the related qq, gg, qg 
fusion subprocesses as the functions of the renormalization/factorization scale Q with nijj- = 175 GeV 
at the Tevatron and mfj- = 250 GeV at the LHC, are shown in Fig.5(a2) and Fig.5(b2) separately. 
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M„.(GeV) M^{GeV) 

Fig. 6. The total NLO QCD corrected cross sections {(Tnlo) and the corresponding relative 
corrections((5) of the processes pp/pp tbH~ + X with niH- = 175 GeV, |r/{,| < 2 at the Teva- 
tron and niH- = 250 GeV, \r]h\ < 2.5 at the LHC, as the functions of mu-. Fig.6(al) and Fig.6(a2) 
are for the process pp tbH" + X at the Tevatron and Fig.6(bl) and Fig.6(b2) for the process 
pp tbH- + X at the LHC. 
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Fig. 7. The total NLO QCD corrected cross sections(crArLo) of the processes pp/pp — tbH" + X at 
the Tcvatron and the LHC, as the functions of the anti-bottom quark transverse momentum cut p\^ 
with niH- = 175 GeV , |%| < 2 at the Tevatron and thh- = 250 GeV , |%| < 2.5 at the LHC. Fig.7'(a) 
is for the process pp tbH" + X at the Tevatron and Fig.7(b) for the process pp — tbH" + X at the 
LHC. 
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Fig. 8. The total NLO QCD corrected cross sections (cttvlo) s-nd the corresponding relative 
corrections((5) of the processes pp/pp tbH~ + X with nifj- = 175 GeV, |r/f,| < 2 at the Teva- 
tron and m^- = 250 GeV, \r]h\ < 2.5 at the LHC, as the functions of tan/3. Fig.8(al) and Fig.8(a2) 
are for the process pp tbH" + X at the Tevatron and Fig.8(bl) and Fig.8(b2) for the process 
pp tbH- +X at the LHC. 
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Fig. 9. The transverse momentum distributions of tlie final state particles(t, b, H~) at LO and NLO 
of the processes pp/pp tbH~ + X with ttlh- = 250 GeV at the LHC and ruH- = 175 GeV at the 
Tevatron. Fig.9(al, bl, cl) are for the transverse momentum distributions at the LHC and Fig.9(a2, 
b2, c2) for the transverse momentum distributions at the Tevatron. 
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